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HOW DID THE IGM BECOME ENRICHED? 

Anthony Aguirre 1 and Joop Schaye 2 



Abstract. The enrichment of the intergalactic medium with heavy ele- 
ments is a process that lies at the nexus of poorly-understood aspects 
of physical cosmology. We review current understanding of the pro- 
cesses that may remove metals from galaxies, the basic predictions of 
these models, the key observational constraints on enrichment, and how 
intergalactic enrichment may be used to test cosmological simulations. 



1 Introduction 

The widespread existence of metals outside of galaxies has been known for a decade 
(e.g., Cowie et al. 1995) by their absorption lines in high-z QSO absorption spectra, 
and was quite surprising when discovered. It has since become clear that this 
enrichment is related to many other aspects of physical cosmology, and that the 
question of "how did the intergalactic medium become enriched?" has a number 
of components that we would like to understand. Among them: 

• The low-density (S = p/ (p) — 1 <J 10) intergalactic medium (IGM), as 
probed by the Lya forest and through Cm, Civ, Sim, Si IV, Ov, Ovi, 
and other tr ansitions, is at least partly enriched at all redshifts and densities 
prob e d (e.g., Songaila and Cowia, 1996t Cowie and Songailal 1998; Elliso n et al 



20001: ISchave et all l2000t ISongailal. l200ll ICarswell et all. l2002t iTelfer et al' 



2002; Bergeron et al., 2002; Schaye et al., 2003; Boksenberg et al., 2003; Pettini et al 



2003 : Aguirre et al.l.l2004lAracil et al.l . l2004ISimcoe et all . 12004 IProchaska et all . 



2004 : iRvan- Weber et all . 120061) . How did these metals get where they are 



? 



• The intracluster medium is highy enriched to Z ~ 0.2 — 0.5^0, with most 
metals produced in clusters res iding in the int racluster medium rather than 

19971 ). How did such efficient metal 



in the galaxies themselves fe.g.. lRen /nn 
ejection occur? 
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Semi-analytic and numerical models of galaxy formation overpredict the stel- 
lar masses of (in particular small) galaxies unless they inco rporate significant 



feedback such as would drive materia l from galaxies fe.g.. lKauffmann et al 



1993; Somerville and Primac kl. ll999f) . What is the relation between enrich- 



ment and galaxy formation? 

The z ~ (and z ~ 2) mass-metallicity re l ation is indicati ve of loss of 
metals by galaxies (e.g., Tremonti et all [20041 : Erb et al. . 200(J ). How much 
baryonic mass and metal mass is ejected by galaxies as they form? 



• Some g alaxies at z 
2(e.g 



Adclberger et al 



(lHeckmanl . l200lh . and nearly all galaxies at z > 
2003), are observed to drive winds that might es- 



cape. Do they escape? What is their impact? 



In connection with these facts, a general picture has emerged that galactic 
winds - driven largely from young and/or starburst galaxies - have enriched the 
IGM both within clusters and the general field. The same feedback may account for 
the dearth of low- luminosity galaxies (relative to the halo mass function) , and also 
the mass-metallicity relationship of galaxies. However, a detailed understanding 
of the various feedback processes is lacking, and there are still open questions, and 
controversies, concerning the time at which various components of enrichment 
occurred, exactly how, and what this tells us about galaxy formation. 

The rest of this article will review recent progress in piecing together this 
picture. First we schematically review the physics behind outflows, the basic ways 
of assessing the predictions of these for the observed intergalactic (IG) enrichment, 
and what general predictions seem to be agreed-upon. We then discuss how the 
metallicity of the low-density IGM is measured, and the basic results obtained 
using these techniques. Finally, we discuss what is probably the best way to get 
a true handle on the enrichment mechanism: direct comparison between detailed 
feedback simulations and the large data sets now becoming available. 



2 Basic physics of galactic winds 



While a single supernova explosion cannot blow material out of any but the small- 
est galaxy, a large set of contemporaneous and co-spatial explosions might. In this 
case, the overlap of super-heated bubbles can overlap and merge into a "superbub- 
ble" that blows out o f the galaxy and forms an outflowing large-scale wind (e.g., 



Heckman et al. , Il990h . A common way to model this has been to assume that 



the superbubble will form a "supershell" for which equations of motion can be 
evolved that include th e driving pressure, gravity, the ram pressure of the medium 
encou ntered, etc. (e.g., Aguirre et all 1200111 iTheuns et al . 2001 : Bertone et al 



20051) . Even if the shell fragments (as is expected), another shell might form, and 
the evolution can be followed by encompassing repeated shell-fragmentation in the 
form of an "entrainment fraction" specifying how much of the ambient material 
must be carried by the outflow along with the superheated gas. 
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A second basic mechanism has more recently been proposed by iMurrav et al 
( 20051) . in which it is momentum depositio n by radiation press ure (coupled via 



200^), that drives the 



dust to the gas3), or cosmic ray pressure ( Socrates et al 
outflow. 

These two possible mechanisms may coexist and would both result in a rather 
complicated outflow structure. A first-order model of the effect of these outflows 
on the host galaxy and the nearby IGM can, however, be generated by assuming 
that whatever the outflow does at very small radii, it eventually organizes into a 
wind with some mass outflow rate m s f r , velocity V w - U1 d, and entrainment fraction, 
and does so in such a way as to conserve some quantity (energy or momentum). 
The two mechanisms differ in the suggested relation between, in particular, m out 
and Kvind- For energy driven winds, it has been argue d that Vwind should b e 



roughly independent of the progenitor galaxy mass (e.g., Heckman et al. . 1990h . 



Then if the feedback efficiency is also fixed, equating the energy generation to the 
energy carried by the wind indicates that: 

TO out (X TO s f r , Krind — COUSt. 

For winds driven by radiation pressure, it is more natural for both the wind 
speed and outflow rate to depend on the progenitor galaxy mass ( or velocity dis- 
persion a), giving a wind possible win d prescription at large radii of (jMurrav et al 
2005HOppenheimer and Pavel [2006h : 



m out oc TOsfr/cr, 1/ W ind — 3er, 

with the wind momentum outflow rh ou tVwind equated to a fixed fraction of the 
momentum input (i.e. luminosity) of the source. 

While these are two interesting perspectives to take, a very great deal happens 
between an erg of energy (or g cm/s of momentum) being deposited, and the 
wind several kpc away. In particular, while energy and momentum are certainly 
conserved in these systems, winds may radiate much of their energy and much of 
the momentum may cancel out if, as in starbursts, it is generated by a distribution 
of sources. 

That being said, suppose we have some prescription for the wind parameters as 
a function of galaxy properties. How might we make predictions for intergalactic 
enrichment? 



3 Modeling galactic winds 

Two basic approaches are used in making predictions for IG enrichment from a 
given wind model: the semi-analy tic, and the numerica l. 



S emi-analytic approaches (e.gjAguirre et all^OOlrilMadau et alll200ll:IScannapieco et al 



Bertone et all 2005 ; Pieri et alT 2006al ) generally proceed as follows: 



2002HThacker et all [2002llFurlanetto and Loebf . 120031 : iBianchi and Ferraral 120051: 



1 If the dust is somehow not coupled to the gas, t h en an outflow of dust coul d ensue and 
contribute to IG enrichment; sec Aguirre ct al. (2001a); Bianchi and Fcrrara (2005). 
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1. A cosmological simulation is used to determine galaxy locations, halo masses, 
and gravitational potentials. 

2. The simulations outputs themselves (for hydrodynamic simulations) or semi- 
analytic prescriptions overlaid on the halos (for dark matter simulations) are 
used to determine galaxy star-formation rates, metallicities, gas properties, 
inflow rates, etc. 

3. A prescription based on a simplified wind model is applied to the galaxy 
properties to determine how far a wind would propagate during a "time 
step" and how the expelled metals would be distributed. These metals are 
then "painted on" to the gas or dark particles. 

4. The metallicity of particles is then tracked, to yield the metallicity distribu- 
tion at later cosmological epochs. 



In the numerical approach (e.g . . | Cen and Ostrikerl 1999 ; Theuns et al. . 2002t 



Thacker et al. , 2002 ; Aguirre et al.ll2005tlCen et all 2005 ; Qppenheimer and Pavel 
20061) . one instead: 



1. Directly simulates gas physics anew for each trial. 

2. Uses prescriptions for star formation, feedback of both energy and metals into 
the gas surrounding star-forming regions, and possibly for the launching of 
winds H 

3. Tracks metals accumulating in the gas particles, and potentially uses these 
metallicities in the cooling rates of the gas. 

Relative to the numerical method, the primary advantage of the semi-analytic 
method is speed: many realizations can be computed on a desktop computer in 
the time that one simulation can be run on a supercomputer. Another advantage 
is that semi-analytic methods allow a clearer connection to be made between the 
physical ingredients of the model and the results, possibly leading to a better 
understanding of the actual physics that is responsible for particular effects. 

Numerical methods, in contrast, treat the gas physics much more accurately 
and can, for example, assess the impact of feedback on the thermal state of the 
gas. Depending on the resolving power of the simulations, fewer prescriptions for 
physics need to be used than in the semi-analytic method, since the physics can be 
directly simulated. Simulations naturally allow a self-consistent treatment of the 
effect of feedback on star formation. Finally, they easily allow for the generation 
of mock observations. 

There is a fairly large literature on IG metals using both approaches, so we 
here focus on what appear to be general and (fairly) well-agreed upon findings: 



2 Of course, in principle, the wind — if driven by energy injection - should launch itself. That 
this generally does not occur without some sort of assistance is generally attributed to insufficient 
resolution. 
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• Galaxies at z > 2 essentially all drive winds: Whether winds are launched 
"automatically" or via a prescription, if the mechanism is consistent with 
winds from low-z starbust galaxies, then it will also predict that most high-z 
galaxies drive them, since the specific SFRs are so much higher then. 

• Wind propagation and escape is quite sensitive to the cntrainment fraction 
and to fwind- This occurs because the two primary forces limiting wind 
propagation are the galaxy's potential well (if the wind speed is not much 
higher than the escape velocity) and the ram pressure of the (potentially 
infalling) gas that must be swept up even if the wind is fast. Moreover, if 
entrainment is significant, then the mass over which the wind energy and 
momentum must be shared may be much greater. 

• The metallicity of the IGM resulting from enrichment by galactic winds is 
highly inhomogeneous and probably has a filling factor <; 10%. 

• Highcr-rcdshift enrichment tends to lead to more homogeneous enrichment 
at lower-redshift. This is somewhat contentious and depends in detail on the 
assumptions made, but appears to be borne out by numerical calculations. 
It occurs because high-z galaxies are relatively smaller and more numerous, 
and because the winds have more time to propagate. It appears to occur in 
spite of the fact that the earliest galaxies of any given mass tend to form in 
highly biased and highly-clustered regions. 

• There is a strong correlation between the gas metallicity and the gas density 
which follows simply from the fact that galaxies form preferentially in high- 
density regions. 

The numerical treatments of IG enrichment reinforce many of these conclusions. 
In addition they seem to generically suggest that: 

• A significant fraction of the metals ejected end up in hot (;> 10 5 K) gas, 
especially when metal-line cooling is not included, but also even when it is. 

• Except at the highest rcdshifts and with the smallest sources, winds tend to 
preferentially propagate outside and in-between filaments, and also have a 
low filling factor, so that there is little overall disruption of the structure or 
statistical description of the Lya forest itself. 

4 Assessing metals spectroscopically 

Nearly everything we know about the enrichment of the IGM at z ^ 2 is via 
optical C iv, C III, Si IV, Si III, O v, and O vi absorption- line spectroscopy of distant 
quasars. Below we briefly review how these assessments are performed, and the 
basic results obtained to date. 
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4.1 Techniques 



The standard technique is to fit line profiles to a set of absorption features to 
derive (where the line is not too saturated) column densities and (where the 
line is resolved) line- widths. This technique works well in the IGM for rela- 
tively high-column lines where there is relatively little contamination. Given 
extremely good data, this can probe gas densities at only a few times the cos- 
mic mean ( Simcoe et al ., 2004). Once determined, the ionic column densities can 
might be converted into metallicities as described below. The prime advantages 
of line-fitting over the pixel-based techniques described next are that it (a) is fa- 
miliar and straightforwardly interpreted, (b) treats high column-density lines well, 
and (c) gives information on line widths. Its main disadvantages are that it (a) is 
unsuitable for statistical searches for absorption that is weak compared with the 
noise or contamination; (b) is subjective in the sense that decompositions in Voigt 
profiles are non-unique; (c) is slow and dependent on prior identification of lines, 
making it unsuitable for application to large datasets and simula t ions. 



T he pixel optical depth (POD) method (ICowie and Songailal 119 98: Songaila 



1998; El lison et al 
2002[lschave et al 



2004 



1999 [|2000l:ISchave et alll2000l:lAguirre et al.l]2002HTelfer et al.l 



. ^OOalAracil et alll2004tlAguirre et alll2004tfpieri and Haehneltl 
Son gailal . 120051 : IPieri et all l2006bl ) is a complementary approach that has 



certain advantages, particularly when applied to weak, widespread absorption that 
may be easily confused with contaminating lines or when applied to large date sets 
or simulations. 

The key steps in the POD search are: 

1 . Compute arrays of apparent pixel optical depths, r app (z) = — In F{z) / -F CO nt {z) 
where F and i^ont are the transmitted flux in the pixel and the local contin 



uum, respectively. Sophi sticated methods have been developed (jAguirre et al 



2. 



2002: ISchave et al.l . l2003l ) to correct for contamination, noise, continuum fit- 
ting errors, etc. in these "apparent" optical depths, taking advantage of the 
fact that many ions generate multiplets (e.g., Hi) or doublets (e.g., ClV, 
Nv, Ovi, and Siiv), with known optical depth ratios. 

Choose a "base" and a "target" transition and bin the pixel pairs according to 
the recovered POD of the former (e.g., C iv as a function of H i). The combi- 



20031 ). and SimCSiiv) (|Aguirre et all l2004h 



natio ns that hav e so far proved to be most useful are C IV (Hi) (ICowie and Songailal . 

1998fl. OvifHi'l (ISchave et al.l.l2000h. SilvfC w) (jAguirre et alll2004h . C ilhC iv) (|Schave et al 



3. Compute the median (or any other percentile) of the target transition POD 
as a function of the binned base transition POD. A correlation reflects a 
detection. The method can be generalized to measure the full distribu- 
tion of target optical dept hs by simultaneously measuring multiple per- 
centiles (jSchave et all 12003) . 



While elemental abundances are desired, we can only measure ion abundances, 
whether from optical depth ratios or column density ratios. This is unfortunate in 
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that uncertainties in the ionization balance are currently the limiting factor in the 
study of intergalactic abundances. But it is also good news, because it means that 
we can constrain the physical conditions in the gas, since the ionization balance 
depends in general on the radiation field, the gas density, and the gas temperature. 

The radiation field (if uniform), can be described in terms of a normalization (in 
the form of the H I ionization rate) and a spectral shape (determined by the type 
of sources, and t he transfer of radiation thr ough the IGM). The UV background 
(UVB) models of Haardt and Madaul ( 200lh have become fairly standard, (which 
is not to say they may not be quite incorrect!) and provide a useful benchmark. 

In photo-ionization equilibrium, the gas density is closely related to the H I 
optical depth: nni oc p 2 T^j, and the Hi ionization rate Thi can be measured. A 
similar (but more detailed) relation can be obtained using mock spectra drawn 
from hydrodynamics simulations. This relation can be used both to perform ion- 
ization corrections, and to indicate the density of gas being probed. 

If the gas is very hot and/or dense, then collisional ionization may dominate 
over photo-ionization, in which case the ionization balance will depend only on the 
temperature. If photo-ionization dominates, then the ion fractions are still weakly 
dependent on the temperature because the recombination rates are. However, as 
long as T ~ 10 4 K, as is usually assumed, the temperature is not the main source 
of uncertainty in the analysis. 



4.2 Observations of IG metals at z > 2 

Using the spectral data and ionization-correction modeling just discussed, in the 
past decade or so we have assembled a reasonably detailed assessment of IG en- 
richment - albeit with some large outstanding questions. 
Basic results thus far include: 



• The carbon abundance (probed via C IV) is inhomogeneous: at fixed overden- 
sity 5 and redshift z, the p.d.f. for [C/H] is a gaussian of width ~0.5- 1 .0 dex , 



i.e. 



the distribution is lognormal (|Schave et al. , 20031: Simcoe et al. , 2004 ). 



Similar results hold for oxygen as probed via Ovi (jSimcoe et all I2004T ) . 



The median carbon a bundance [C/H] (an d also the width of the p.d.f.) in 



creases with density (jSchave et all 120031 ). and there is (so me) carbon in 
som e) underdense gas (corresponding to Ahi ^ 3 x 10 13 ) (jSchave et al 



2003) 



For a standard choice of the UV backgroun d (HM01 Q+G UVB ) , the median 
[C/H] does not evolv e from z ~ 4 to z ~ 2 dSchave et al. . 120031). In addition 



meta l s exist at z >, 5 (|Songailall2001l : lPettini et al 
l2006l : ISimco3 . l2006l) . 



2003 



iRvan- Weber et al 



• The I GM is alpha-en hanced relative to the sun: [Si/ CI 
20041 ). [O/C] £ 0.0 (|Telfer et all 120021 : ISimcoe et al. 



' 0.2- 1.5 (jAguirre et al 
2004h . This suggests 
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enrichment by predominantly Type II supernovae, perhaps even with a con- 
tribution by hypernovae (see, e.g., Qian and Wasserburd . 2005 . and other 
contributions to this volume). 



Results are UVB-dependent. A lHaardt and Madaul (|200lh UVB model with 
galaxies and quasar light included, and 10% escape fraction from galaxies 



seems plausible; a much harder or softer does UVB does not (jSchave et al 
120031 : lAguirre et all l2004t ISimcoe et all 12004 ) 



Using detected Sim and Cm absorption, it can be sho wn that most strong 



Si IV and Civ abso rption arises in photoionized gas ( Schave et al . 20031 



Aguirre et all 120021 ) 



• Metals are correlated with themselves and with g alaxies (jAdelberger et al 
20051 ; iPieri et all l2006tj ; IScannapieco et all |200< 



4.3 Basic Implications of derived abundances 

A basic implication of the above findings is that essentially all of the results are 
explainable by metals "sprinkled" onto simulations with no winds, at z p> 5. But 
does it mean that all metals in the IGM come from such an early enrichment 
period? There are some reasons to think not. 

First, a substantial fraction of all metals generated by stars at z > 3 appear to 

reside in the diffuse IGM: [C/H] 2.8, [Si/H] ~ - 2.0, [O/H] 2^3 (for HM01 

Q+G) imply /zjgm 10 - 30% (jBouche et all |2006[) . 

Second, because their ionization fractions fall quickly with temperature above 
~ 10 5 K, Civ and Siiv cannot probe very hot gas - metals in such gas would 
simply be invisible to probes using these ions. 

In general, the rough amount of metals present in the IGM can be reproduced 
by models that include galactic winds (though not with models that do not include 
winds). However, the observations summarized above present a detailed and chal- 
lenging target for the next generation of models. In the next section we summarize 
some recent work comparing cosmological simulations - including enrichment - in 
detail to the observed spectra. 



5 Confrontation with simulations 

While much can be inferred using observations directly, It is very useful to compare 
them to synthetic spectra generated from simulations with enrichme nt - including 
the eff ect of the outflows on the IGM. In one investigation of this sort . lAguirre et al 



( 2005f ) analyzed two sets of h i gh-resolution particle-pase d hydrodynamic simula- 
tions (jSpringel and Hernquistl 20031 Theuns et al. . 2002f ) with different feedback 
prescriptions that both drive strong winds at z J; 2 as per the late enrichment 
scenario. Four general and qualitative conclusions of this study were as follows. 

First, the feedback simulations produce far too small values of tciv/thi and 
tciii/tciVj f° r & H UVB models employed. 
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Second, metal-rich gas in the simulations is too hot and too low-density: es- 
sentially all of the intergalactic and enriched gas in both simulations exists in 
low-density, high-temperature (10 5 K <, T <, 10 6 K) bubbles. Since both the 
Civ/C ratio and the Ciii/Civ ratio fall off quickly with both increasing temper- 
ature T <; 10 5 K and decreasing density 5 ^ 10, the gas becomes nearly invisible 
in C iv, and where C iv is detected there is virtually no accompanying C III. 

Third, metal line cooling, which was not included, should be important. If 
crudely modeled (by assuming that gas cools to T « 10 4 K if its cooling time 
is shorter than the Hubble time), much more Civ absorption is present and the 
median tciv / Thi could plausibly match the observed one. 

Finally, metals nonetheless appear too low-density and too inhomogeneous: 
While including metal cooling might lead to a sufficient increase in C IV absorption, 
it does not appear to fix the difficulty in reproducing C ili/C IV. This can be traced 
to the fact that the metal-rich gas has too low density to produce the observed 
ratios. This problem may also be lessened by a proper treatment of metal cooling, 
as the cooling would also affect the gas dynamics, allowing gas to contract as 
it cools. But there is yet another problem, which is that although the median 
tciv /thi can be roughly reproduced by the simulations, the spread in tqw at 
a given thi cannot: the metals in the simulations are too inhomogeneous. This 
problem seems unlikely to i mprove from a correct tr e atmen t of metal cooling. 

In a more recent study, lOppenheimer and Pavel (|2006f) have compared simu- 



lated spectra from a large suite of feedback simulations to observed line and pixel 
data. (See also contributions by Dave and Oppcnheimer in the present volume for 
more details; here we just mention a few points.) 

The Oppenheimer & Dave simulations include metal-line cooling, and wind 
prescriptions chosen to correspond to "momentu m driven" and "energy driven" 
winds as discussed in Sec. O In comparison to Aguirre et al. ( 20051 ) there are 



two findings of note. First, Dave et al. find much better agreement in Civ/Hi 
an d Chi/Civ op t ical d epth ratios between their simulations and the data than 



do lAguirre et al.l (|2005f ) . This is very likely due to the inclusion of metal line 
cooling. Second, in models for which the C iv optical depths are compatible with 
ob servations, th e over all cosmic carbon abundance is higher than that inferred 



bv ISchave et al. (2003) based on the same data. This indicates that (as argued 



by Dave et al.) metals are being hidden in relatively hot gas not probed by the 
observations. 

It would be very interesting to see how the Dave et al. simulations fare when 
applied to O vi, which is much more sensitive to the presence of hot gas, and when 
applied to percentiles in C iv other than the median - which Aguirre et al. ( 2005h 



found difficult to reproduce using simulations, even with cooling included. 
6 Concluding facts and opinions 

We will conclude with some facts, and some opinions. The key facts are: 

• Both observations and simulations of IG enrichment are getting drastically 
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better, and the former are, finally, strongly constraining the latter. 

• Comparisons between observations and simulations have to be done with 
care. 

• Observations of IG metals strongly constrain models of feedback from star 
formation (a crucial ingredient of simulations and galaxy formation models). 

Some opinions: 

• Enrichment is unlikely to be purely "contemporary" or "primordial" : En- 
richment almost surely occurs at z < 4, perhaps mostly hidden in hot gas. 
However, it is unclear that z < 6 enrichment can account for all measure- 
ments. 

• It seems likely that a galaxy-mass-dependent V w i n d is necessary - but it is 
not clear that this requires radiation-pressure-driven winds. 

• There is a dearth of high-quality, high-resolution simulations of individual 
wind-driving galaxies in a cosmological background. Please do some! 
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